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Abstract 

Kontsevitch's work on Airy matrix integrals has led to explicit results for 
the intersection numbers of the moduli space of curves. In a subsequent 
work Okounkov rederived these results from the edge behavior of a Gaussian 
matrix integral. In our work we consider the correlation functions of vertices 
in a Gaussian random matrix theory , with an external matrix source. We 
deal with operator products of the form < l\i=i jjtiM ki >, in a 4 expansion. 
For large values of the powers hi, in an appropriate scaling limit relating 
large fc's to large N, universal scaling functions are derived. Furthermore 
we show that the replica method applied to characteristic polynomials of the 
random matrices, together with a duality exchanging N and the number of 
points, allows one to recover Kontsevich's results on the intersection numbers, 
through a simple saddle-point analysis. 
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1 Introduction 



Random matrix theory (RMT) has been applied to many physical prob- 
lems, and also to mathematical subjects such as the distribution of zeros of 
Riemann zeta function or combinatorial problems and it has led to several 
meaningful results [1] . It also plays an essential role in the theory of random 
surfaces and for string theory. Several kinds of correlation functions in ran- 
dom matrix theory have been studied. In previous papers, we have studied 
the correlation function of the eigenvalues [4], and the correlations of the 
characteristic polynomials [7, 9], for which we have derived explicit integral 
representations. 

In this article, we consider the correlation functions of vertices on the 
basis of previously derived integral representations. The diagrammatic rep- 
resentation of the vertex (trM fc ), where M is a random matrix, is obtained 
through Wick's theorem, by the pairings of fc-legs, each leg carrying the two 
indices of the matrix element M^. For N x N matrices, the two indices 
run from 1 to N : i,j — 1, N. 

We restrict ourselves in this article to complex Hermitian random matri- 
ces. The distribution function for M is Gaussian with an external matrix 
source A. 

P A ( M ) = — e-T ^-NtrMA (1) 

Za 

When one sets A = 0, it reduces to the usual Gaussian unitary ensemble 
(GUE). 

The correlation functions for the vertices V(ki, k n ) are defined as 

V(ki, k n ) = < trM kl trM k ' 2 ■ ■ -trM fc " > (2) 

The normalization is chosen so that they have a finite large- N limit. These 
functions are closely related to the Fourier transform of the correlation func- 
tions of the eigenvalues, 

/•oo „ N 

U(t 1 ,...,t n ) = / e t ^ ux >R n (\ u ...\ n )Y[d\ l (3) 

J -co y 

where the correlation function of the eigenvalues is 

RniX,, A n ) =< I] ^tr5(A, — M) > (4) 
i=i iV 

Indeed 

1 n 

U(t 1 ,...,t n )=<—Y[tre^ M > (5) 
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are generating functions of the V{k) since 



u{t u ...,t n ) = < trM ' 



ki=0 



oc 



kl trM k2 ■ ■ ■ 



trM kn > 



{tti) kl ■ ■ ■ {it n ) kn 
k 1 \k 2 \---k n \N n 



(6) 



When the distribution of the random matrix is Gaussian, the average 
of the vertices gives the numbers of pairwise gluing of the legs of the vertex 
operators. The dual cells of these vertices are polygons, whose edges are pair- 
wise glued. We thereby generate orientable surfaces, which are discretized 
Riemann surfaces of given genus. 

Okounkov and Pandharipande [16, 17] have shown that the intersection 
numbers, computed by Kontsevich [15], may be obtained by taking a simul- 
taneous large N and large ki limit. Furthermore the correlation functions of 
these vertices are interesting, since they give universal numbers in the large 



We have investigated in an earlier work the F.T. of the n-point correla- 
tion function U(t 1 ,...t n ) for the GUE, and found a simple contour integral 
representation valid even for finite N [3, 4]. 

In this article, we extend this integral representation to the vertex corre- 
lations V(ki, ..k n ), and examine the scaling region for large ki and large N. 
In this integral representation, the asymptotic evaluation by the saddle-point 
method requires a careful examination to deal with pole terms. This leads 
to a practical way to compute intersection numbers which we discuss in de- 
tail. We also show that the F.T. of the correlation functions (C.F.) of GUE 
near the edge point of the support of the asymptotic spectrum, is equiva- 
lent to Kontsevich's Airy matrix model ; the identification is based on the 
replica method and over a duality for computing averages of characteristic 
polynomials. 

The article is organized as follows. 

In section 2, we consider the F.T. of the one point correlation function at 
a bulk generic point in the large N limit. This is done by a contour integral 
representation, and we obtain the behavior of < trM 2k > when N and k are 
large. We show that in this limit, one recovers the behavior of the one point 
function near the edge point of the spectrum. 

In section 3, we consider the correlation function of two vertices. 

In section 4, we investigate the correlations of the n-vertices. 

In section 5, we introduce a replica method, relying on averages of char- 
acteristic polynomials. This, together with a duality, allows us to make con- 
nexion with the Kontsevich model, recovering thereby generating functions 
for the intersection numbers. 

In section 6, we present a short summary. 



N limit. 
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2 One point correlation function 



The correlation function i? n (Ai, A n ) defined by 

R n (X ± , X n ) =< f[ ^tr(5(A, -M))>. (7) 

is thus equal to 

Rn{K A n ) =< ft ^tr( ^e-*^-*) > . (8) 
The Fourier transform U(ti, t n ) of R n (Xi, A n ) is thus given by 

U(h,...,t n ) = / HdX i e i ^ tiXi R n (X 1 ,...,X n ) (9) 

J-oo i=l 
i n 

= <^n^ M >. (10) 

Note that U(ti,...,t n ) is normalized to one when all U = 0. The function 
U(ti, t n ) is the generating function of the correlation V(ti, t n ) as shown 
in (6). 

These F.T. of the correlation functions U(ti, t n ) were investigated in 
our earlier study of the kernels for the correlation functions [2, 3] ; there we 
had derived an exact integral representation for these correlation functions. 
We consider the probability distribution 

P A ( M ) = — trMi-NtrMA 

Za 

then one finds the exact result 

U( tl ,...,t n ) = —±— e -&tt+"^ 

(ti ■ ■ - tn) 



X 



n N it 



(2ni) n A = AA = AV N(u- ai ) 



x i U P M 9 + AT fa tq)]( U p U q) 

where the integration contours circle around the eigenvalues a 7 , (7 = 1, • • • JV), 
of the source matrix A in the anticlockwise direction. 
For the one-point function 

U(t) =< -^tr e ltM > (12) 



the exact integral representation for finite N is thus [2] 

which reduces, for the pure Gaussian model a 7 = 0, to 

U(t) = - I + -iL) N e-^ +itu . (14) 

w it J 2m y Nu v ' 

For this sourceless GUE one obtains immediately in the large N limit 
The generating function for Bessel functions Jj(x), 

oo 

e Hu + l) = £ {m y M 2t) 

3=-oo 

J_,-(x) = (-l)Vi(x), (16) 

leads to 

tf(t) = iji(2t). (17) 
The semi-circle law for the density of states of the GUE follows : 

dt 



/oo fit 



= o, (£>i) 

i^(f?, (f < I)- («) 
Returning now to the exact expression (14) for finite N one finds 



f Nk L ^ r ( Ar -^ + 0r(A;-/ + i)r(A;-/ + 2)r(/ + i)2' J 



(19) 

There are no odd powers of in this expansion, as is well-known for the GUE 
case, for which the successive terms of the 1/iV-expansion are of the form 
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1/N 29 , where g is the genus of the surface generated by the Wick contractions. 
From the relation between U(t) and < trM 2fc >, we obtain 

1 (2*;)!* T(N) 

-<trM 2fc >=^f[£ 



N ' N k l f^T{N -k + l)T{k-l + l)T{k-l + 2)T{l + l)2 li 

(2k)\ (N — 1)(N — 2) • • • (AT — fc) (N - 1)(N - 2) ■ ■ ■ (N - k + 1) 
iV fe ^ fc!(fc + 1)! + (fc- l)!fc!2 

(jV-l)(jV-2)-(jV-fc + 2) 
+ (ifc-2)!(Jb-l)!8 + 1 [ ' 

This exact representation leads to (20) the expansion 

I < trM2 , (2ib)! fc(fc-l)(fc + l) 

iy <trM >-^ !( ^ + 1) ,L i + 12iy2 

fc(fc + l)(fc - l)(fc - 2)(fc - 3)(5fc - 2) 1 
+ 1440AT4 + ^6^ ^ 

The large N limit is the first term of the above expansion : 

l im _L < t rM 2k >= ,.j\ k) \.. (22) 
n^ooN fc!(ife + l)! 

the k-th Catalan numbers. For large k this number behaves as A=\A k . 
Therefore the resolvent 

has a square root singularity at \ 2 C = 4. This corresponds of course to the 
vanishing of the asymptotic density of state (18) as a square root at the edge. 
Returning to the expansion of < trM 2k > one finds for large k 

k{k 2 -l) fc(fc + l)(fc - l)(fc - 2)(fc - 3)(5fc - 2) 1 
[i+ 12iV 2 + 1440iV 4 +U{ N* )V 

(24) 

2 

When k is of the order k ~ iVs, the above series exhibits a scaling behavior: 

1 1 U3 L6 I 

iV <trM2t> ~^ 4l|1+ l^ + (12P^ + M' (25 » 
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The power 4 fc corresponds to the location of the edge of the support of the 
asymptotic spectrum (A 2 = 4), and it is not universal. But the successive 
terms being powers of k 3 /N 2 is universal, since this feature is related to 
the square root vanishing of the density of states. The scaling function in 
this double limit of large N and large k when k behaves as k ~ iV 2//3 , is 
also universal. The above coefficients of the scaling function in (25) provide 
the intersection numbers of the moduli of curves. The universality of the 
coefficients of the series in powers of k 3 /N 2 corresponds to the F.T. near the 
edge of the spectrum A ~ A c , for the universal Airy kernel. 

In order to obtain the full scaling function of k 3 /N 2 , and not simply the 
first terms of the expansion as in (25) , we now consider an exact integral 
representation for < trM 2fc >. From (6) 

1 < trM 2k > = (2*)! 1 dt m 



N v '' J 2nii 2k t 2k+1 



(2k)\ r r dtdu 1 it 



1 + —)N e ~ m +itu ( 26 ) 



i{-l) k J J {2ni) 2 t 2k+2K Nu 

For large k, large N, we apply the saddle-point method to (26). The 
integrand in (26) behaves as e N< ^, with 

it 

<p = -2fclnt+ - + itu (27) 
u 



it[-— + 1] = 



The saddle point equations, 

1 

du 1 u 
d<f> 2k .A . n 

give as solutions u = ±l,t = ^fik. Expanding around the saddle-point in a 
standard way one obtains in the large k, large N limit 
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1 . k 2k !e 2k 1 k3 4 k 

— < trM 2k >~ \ ' : : ^e^ 29 

N k 2k +! 4?rk 0fkl 

We have thus obtained the scaling function 

k 3 k 3 
/(— ) = exp[^— ] (30) 

in accordance with the expansion found hereabove in (25). 

Instead of the bulk spectrum, we now consider the F.T. of the one point 
correlation function near the edge point (A ~ A c ). We denote by U(t) the 
F.T. of the one point correlation function near the edge of the spectrum, to 
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distinguish it from the bulk U(t). We expect from the previous argument 
that this U (t) becomes < trM 2k > if one puts t = —ik. 

To explore the vicinity of the edge, it is convenient to introduce a trivial 
external matrix source whose eigenvalues are all a 7 = —1 in (13), and mul- 
tiply e~ lt in order to compensate for this uniform shift. Thereby the edge is 
now at the origin and 

e lt f du it s N *L +itu 



U(t) = + — -)"e~™ 

w it J 2-ki N(l+u^ 



-■it 



du e M«a+tffe) e -&+«« (3i) 



it J 2ni 

In the regime in which t ~ N 2 / 3 and u ~ N~ 1 / 3 one may expand for small 
u up to order u 2 , and the contour integral becomes a saddle point Gaussian 
integral. Note that the term itu is cancelled in the exponent. We have 

T-T^S 1 mt [°° (III itu 2_± t 2 u 

U(t) = —e^ 7 / e ^ 

it J-oo 2iri 

1 (it) 3 

el^FT (32) 



We thus recover the scaling function / = exp[k 3 /12N 2 ], when we put t = —ik 
in (32). The difference between -^ < trM 2fc > and U(t) when we replace it by 
k, is only the prefactor |4 fc , . We have thus shown that the scaling function 
of (30) in < trilf 2fc > for k ~ iV 3, namely exp[ 1 ^ v - 2 ], may also be obtained 
from the F.T. of the one point correlation function near the edge by setting 
t = —ik. 

Changing the prefactor by multiplying by l/t 1//2 , and defining x = —it/2 1 / 3 , 
we obtain the generating function, 

1 3,3 
Fix) = — e24Af2 

X 1 

1 X x 4 

x 2 + 24N 2 + (24) 2 2!iV 4 + " ' 

= E<^- 2 >^ 3 - 2 . (33) 

The numbers < > in this expansion coincide with the intersection 
number of the moduli of curves, as we will be justified below by the replica 
method. 

From the expansion (33), we obtain < Tj > as 

<^>*=(24^i (2 = 0,1,2,...). (34) 
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These numbers agree with the values of the intersection numbers computed 
earlier by Kontsevitch and Witten [15, 19]. For < r >, we need a special 
consideration in (33). If we put g — |, we get < r >. However by definition 
< jjtiM 2k >= 1 for k — 0. Therefore we define < t > 9 =o= 1 instead 
of < r_ 2 > 9 =o= 1 for comparison with the intersection numbers. We will 
discuss these intersection numbers later. 

We have used here the integral representation to derive U(t). Since this 
is related to the edge problem, we could have used instead the Airy kernel 
K A (X,X), 

K A (\,») = ^WAfc) ~ A(A)4j 

Ai(X + z)Ai(n + z)dz (35) 



r 

Jo 
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where the Airy function A^x) is given by 



Mx) = — / ( ^ -"v/C (36) 

Z7T J-oo 

Near the edge of the spectrum A = 2, the density of state is given by Ka{\ A) 
in the appropriate large N scaling. Let us verify that one can recover the 
previous result from there : 

U(t) = / d\e itx / dzA^X + z) 

J-oo JO 

I OO P OO P OO /^OO 

= " dz e i*A A e K 3 +^(A+z)^ /• e ^-H(- t+ *)„ d 

(27r) 2 JO J-oo J-oo J-oo 



/*oo /* oo /*oo 

JO J— oo J— oo 



| /"OO /'OO /■oo 

2n Jo 

r°o 1 (it) 3 . 
/ — 7=e 12 i2t 
Jo 2Vi7rt 



1 1 



e^ (i ' )3 (37) 



2v^F (ii)3/ 2 
which coincides as expected with (32). 

3 Two-point correlation function 

In the case of the two-point correlation function, we have 

i*2(Ai, A 2 ) =< ^tr<J(Ai - M)^tr5(A 2 — M) > (38) 
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and the F.T. of R 2 (Xi, A2) is 

U{hM) =< ltre^ M ltre^ M > (39) 

This correlation function has been obtained in closed form, for finite N, with 
the help of the HarishChandra-Itzykson-Zuber integral, in [3] 

T j, s v-^ llj<j[ a i ~~ a j + Jftl{^i,ai ~ 3j,ai) + ^2(^1,02 ~~ ^,0:2)] 

U{t 1 ,t 2 ) - 2^ ff ( n _ n \ 

ai,a 2 lli<j\ a i a j ) 

x e itia ai +it2a c , 2 -^tl-^tZ-±t 1 t 2 5 ai , a2 ^q-j 

This sum is then divided in two parts; a.\ = a 2 and a.\ ^ a 2 . The first part 
gives 

U^ti t 2 ) = XI II ~ ~ ^ + + t2 ^ Si,ai - Sj ' ai ^ c~^ (tl+t2)2+iitl+t2)aa ' 

ai i<j ( a i ~~ a j) 

(41) 

This may be expressed as the contour integral 

U\ tl M) = , 1 r I P- [1 + lu 1 + t2 \y ^ u+ ^^ 2 +^ (42) 
V ; i(fi + 1 2 ) 7 2m £J l N(u - a 7 ) J v ; 

which is nothing but U(t\ + £2) from (13). 

The second term is expressed by the double contour integral 

U n (U h) = e -^ t2+ ® l &l(ii ' 2 e"'»' +i '^ TT(1 + — )(1 + — ) 

[tut2) f (2™) 2 T 1 = V i + 7V( Ul -a 7 ) ni+ 7V( U2 -a 7 ) J 

x 1 (ui-u 2 + ^(ih-it^jm-u^ 

hh {u 1 -u 2 + jft 1 ){u 1 -u 2 -jft 2 ) 



Noting that 

1 - 



t\t 2 _ (ui -u 2 + jj(it! - it 2 ))(ui - u 2 ) 



iV 2 (wi - u 2 + jjh)(ui -u 2 - jjit 2 ) (ui - u 2 + jjh){ui - u 2 - j^t 2 ) 

(44) 

we find U H is a sum of a disconnected part and a connected part. 

Therefore, we find the connected part of U n (ti, t 2 ) as the contour integral 



X _|_ l l± ) N (l + ^ 2 ^ N 



(2m) 2 v iVui v Nu 2 

N\ Ul - u 2 + f ) ( Ml - M2 -f) 



(45) 
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We have set here all the a 7 = to deal with the pure GUE. The contour is 
around u\,u 2 = 0. This U c (ti,t 2 ) may be expanded in powers of t\ and t 2 . 
Together with the exponential factor, it yields the expansion 

U^(t u t 2 ) = 1[1 - + t 2 f + itita + 1(1 + -L)(f« + 4) 

+(\ - ± + + tl tl) + ^)t\t\ + O(0] 

(46) 

One may wonder why we have obtained odd powers of in this expression 
but one can check that combined with U{t\ + t 2 ), the odd power in of 
U^ I (t\,t2) cancels. 

If we compute the integral representation of U^ I (ti,t 2 ) by deforming the 
contours to collect the contributions of the poles at Ui = 00 and w 2 = 0, 
instead of u\ = 0, u 2 = 0, we obtain 

U c (h, t 2 ) ( oo,o) = l[*i*2 - \t\t\ - \{t\t 2 + htl) + 0(f)] (47) 
From the expression of U(t) and U* I (ti,t 2 ) found in (19) and (46), we have 

ltf(fi + f 2 ) - t£'(fi, t 2 ) = -1[M2 - 5(^2 + - ^ + 0(t 6 )] (48) 

which is indeed equal to — ^2)(oo,o) m (47) and justifies the deformation 
of the contour integration in order to collect the residue of the pole at u± — 00. 

Let us first consider the connected part of the two-point correlation func- 
tion in the large N limit. 

Using the Bessel function formula (16), we find from the residues at U\ — 00 
and u 2 = 0, 

1 00 

f/ c (ti,t 2 ) = ^E(- 1 ) / ( / + 1 )^+i(2ti)Jm(2t 2 ) (50) 
Noting that the F.T. of the Bessel functions Jj+i(2£) is 

/oo 1 
dtJi +1 (2t)e ltx = 2(4- A 2 )-5 C os[(Z + l)0] ((Z + 1) is even) 
-00 

/oo , 
dtJ l+l {2t)e ltx = 2i(4- A 2 )~2 s in[(Z + 1)0] ((Z + 1) is odd) 
-00 

sin0 = ^. (51) 
10 



we obtain by inverse F.T. 

1 r 00 
ik(Ai,A 2 ) = N J dt l( it 2 e Mi ^ £(-!)'(/ + lM + i(2tiM + i(2t 2 ) 



1=0 



N 2 {2mf 
1 00 

]V5(2^2 E(-l)'(' + i)^* 1 + (-l) m e-^] 

x j e i(j+i)02 _|_ ^^y+ig-^+i)^ 1 ! 



(4-A 2 )(4-A|) 
1 1 4-AiA 2 

(Ai - A 2 ) 2 ^(4 - Af ) (4 - A 2 ,)' 



(52) 



This result agrees with the earlier derivation of [10]. 

We now consider the integral representation for < trM fcl trM fc2 >. We 
have 

-L < t,M*t,M* > c = g^gMl / / ^|||^ (53) 
N 2 (_i)*i+fca J J (2m) 2 tl kl+1 tl k2+1 

where f/ c (ti,t 2 ) = £# 7 (*i,t 2 ) ~ U(t 1 + t 2 ). 

< trM 2fcl trM 2fc2 >^ (2/ci)!(2fc 2 )! / ^1^2 / d«idw 2 e~ 



^2 " 7v 2 (-l) fe i+ fc 2 / (27ri) 2 7 (27ri) 2 tf i +1 t 2fc2+1 

(«l-«2 + f)(«l-«2-f) ^ ^ 

We are interested in the large N and large k\,k 2 behavior, but in the 
region in which the ki are of order iV 2//3 . As for the previous calculation of 
jj < tiM 2k >, after exponentiation , we find again that the saddle points 
are t lc = —ik\ and t 2c = — ik 2 with u± — u 2 — 1 in (54). Then we expand t 1 
and t 2 near the saddle-points 

ti = -i/ci(l + vi) 

t 2 = —ik 2 (l + v 2 ) (55) 

and expand for Vi,v 2 small. The integrations over Vi,v 2 become Gaussian, 
and they are equivalent to the replacement of ti by their saddle point values 
—iki (i=l,2). Therefore we have for large k and N, after the shift Ui — > l + Ui, 

At < trM 2fel trM 2fe2 >" 
N 2 

„ n , s 1 (fe 3 +fc 3 ) f°° d Ul du 2 e kiui + ^u 1+ k 2 ul + ^u 2 

= -C(ki,k 2 )e3N^ ( 1 2 ' / . r- , 7— rr(56) 
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where the constant C(k±, k 2 ) is 

(2k 1 )\(2k 2 )\e 2k ^ 2k - 



C(k 1: k 2 ) 



(-l) k ^^k 2 kl kl kf 2 

4 2fc i+ 2fc 24 7 rA; 1 A; 2 ^ 



The integration in (56) requires a careful examination of the pole terms. 
Since the denominator may vanish, we use 

= P( l — tet) + ™S(ui - u 2 + h (58) 



Mi - u 2 + ^ - «e u 1 -u 2 + % N 



The 5 function contibution is nothing but < trM 2kl+2k2 >. The principal 
part is evaluated by writing the denominator as 

= -i / dae i{ui ~ U2+ ^ )a (Imh > 0) (59) 

u 1 -u 2 + ^ Jo \ i j v ) 

The integration in (56) becomes 

z 2= r duidu2 / oo rfQ;rf/ 3 e EL 1 (^«f+l^)+^(«i-« 2 +^)H-^ 1 -^-^) 

j-oo (27ri) 2 7o 

1 r°° u,( fc l I j(g±gh 2 W fc 2 '("+/3) n2 i tofel »/3*2 „ 

In this representation we have assumed that Imk\ > 0, Imk 2 < 0, but after 
integration over a,/3, this condition becomes irrelevant. We replace 

a + (3 = x, f3 = z (61) 

with 

< z < x (62) 

Thus we write I 2 as 



e an 1 f°° f x 

, / dx I dze 

^ix\Jk\k 2 Jo Jo 



kl+ k 2 . r 2 | i(k 1 +k 2 ) i(k 1 +k 2 ) , 
4,k 1 k 2 2iV iV 



+ 



47r v / M^(A;i + k 



2j 



ik -\ k- 



1 K 2 

— fc l+ fc 2 2 

x / ifcifc2 rfa:e 4fc i fe 2 (63) 
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After mutiplication by e ; W e obtain 



e3jv 



vrCfci+fe,)^ = _ V 2 £ (fcifc 2 (fci + fc 2 ))' 



27r(fci + fc 2 ) ^l!(2l + l)4'iV 2i 

(64) 

Finally 

(65) 

Note that in this evaluation, we have already subtracted the 5-function term, 
when we have represented the pole terms in terms of integrals over a and j3. 

For making contact with Kontsevitch normalization [15] , we change k\, k 2 
to ^J73, 7^73, and multiply a factor J^- . We obtain the intersection numbers 
for two points (n=2) as an expansion of the error function. 

tp, \ 1 (xi+x f " [x 1 x 2 (x 1 + x 2 )] l (-l) 1 

where we set Xi = kj within the appropriate factors. This function F(xi,X2) 
is a generating function of the intersection numbers. If we expand it in powers 
of x\ and x 2 

F(x u x 2 ) = Y,< r h T h > g -£«f (67) 



X\ x 2 



the coefficient < ti x t\ 2 > is the intersection number for genus g ; the genus 
is specified by 3g — 1 = h + l 2 , and then the coefficient is then given by 
returning to (66). 

As seen previously for < trM 2fe >, the asymptotic evaluation of 
< trM 2fcl trM 2fc2 > is given by the Fourier transform U(ti,t 2 ) near the end 
point. 

From the above expression we obtain coefficients, which are the intersec- 
tion numbers for n=2, 

1 1 

< T 2 T > g=1 = < T T 2 > g =l— — , < r 5 r > g=2 = < T T 5 >g= 2 



24 , -.o<u-s=z -'U'0-g=z ( 2 4) 2 2!' 

2 1 1 29 

<r 1 > g =i=u> <r.r 1 > g ^=— 1 <t 3 t 2 > s = 2 = — (68) 

These numbers agree with Witten's earlier results [19]. For < Tq > g =o, we 
use the normalization, < Tq >= 1, in analogy with < r >= 1. 

It is worth noticing that when we set x 2 = 0, i.e. k 2 = in (65), we do 
obtain 

1 x ' A 

F(x l7 0) = F( Xl ) = — e^, < r h r >=< r h > . (69) 

X\ 

The above relation is the string equation, as will be explained later. 
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4 The n-point correlations 



The Fourier transform U(ti, ■■■■,t n ) is given by (11). Using Cauchy determi- 
nant formula, it is expressed as a determinant. For the connected part, we 
take the longest cyclic rings for the indices p, q. 



u c ( tl ,..., tn ) = i 7re -A«;+-«5) 



x TT -i , (70) 

cycle N P P a 1 

where the last product is the maximal cycle for the indices (p,q). For in- 
stance, in the case n=3, we have two longest cycles (1 — > 2 — > 3 — > 1) and 
(1_> 3 -2-1), 

n * = * 

cycle ¥p + Up~U q ( Ul - U 2 + f )(« 2 " «3 + f ) (u 3 ~ «1 + f ) 

+ K - «3 + f )(^2 " U X + f )(% - U 2 + f ) (71) 

These two terms contribute to the connected part of U(ti,t 2 ,t 3 ) as 
U e (t u t2,t 3 ) = ^e-™ { * 1+ *+® 



x 



Nu 3 ' 



r du ld u 2 du 3 clTt ^ (1 + ^(1 + ^(1 + ^ » X 

J (2«) 3 V - U 2 + f )(« 2 - U 3 + f )(« 3 - Ul + f ) 

/i | iti \N(-\ I ita \ ATZ-i i if 3 \N 

(«! - « 3 + f )(«2 - «i + f )(«s - «2 + f Y 



We find that the two terms are identical, and they are symmetric polynomials 
of the t p . 

As we have seen, we have to add several terms to obtain the connected 
part of U(ti, t n ), to deal with the cases ctj = aj in the summation implied 
by (40). The evaluation of this integral is an extenstion of the previous study 
of I 2 . 

For the n=3 case, we have in the large N limit, neglecting all terms, 



du\du 2 du 3 e z 



(27T2) 3 (Ui - U 2 )(U 2 - U 3 ) (U 3 - Mi) 
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du i 1 V- ^yi^w s /W3y 2 

(27ri) 3 ^ 2 J^ 3 V V V 

oo 

x X! ( m i) Sl (»2) S2 (»3) S3 ^ (2ti) J S2 (2t 2 )J Sa (2t 3 ) 

S1,S2,S3 = — oo 

/ ^\mi+m2+m3^2mi+Zi+i3+1^2m2+ii— i2^2m3+i2+i3+l 



E 



m^mj + 1 + Zi + Z 3 )!m 2 !(m 2 + Zi - Z 2 )!m 3 !(m 3 + Z 2 + Z 3 + 1)! 

(73) 



We consider the coefficients of (-if^+HfHfHf 3 of J, which is denoted 
by / 2fc i'2fc2,2fc 3 . We put 

, I/1/3 
m i — k\ — - — — — — 
2 2 2 

^ l ^ J 2 
= fc s ------ (74) 

If Zi is even, then Z 2 is even and Z 3 is odd. If Zi is odd, then Z 2 is odd and 
Z 3 is even. These two cases give the same result, and give a factor 2 for Z x 
even. We change Zi — > 2Zi, Z 2 — > 2/ 2 and Z 3 — > 2Z 3 + 1. 



j2fci,2fc 2 ,2fc 3 

00 1 

= y 

h Ms=o - h - l * - + h + k+ - Zi + Z 2 )!(£; 2 + Zi - Z 2 )! 

X (A; 3 -/ 2 -/ 3 -l)!(A; 3 + / 2 + / 3 + l)! (75) 
This sum is expressed by the contour integration, 

j2fci,2fc 2 ,2fc 3 

_ 1 r dxdydz (l + x) 2kl (l+y) 2k2 (l + z) 2k3 

" (2fci)!(2fc 2 )!(2fc 3 )! * (2yri) 3 x^+^l - X y){l - xz){l - f) 
1 rdxdydz (x + y) 2kl (1 + y) 2k2 (1 + zy) 2ka 

" (2fci)!(2fc 2 )!(2fc 3 )! J (2vri) 3 ^1^1+^3+1^3(1 - x)(l - *)(1 - xzj 

(76) 

where the contours are around x — y — z — 0, and in the last line, we 
have made the change of variables x — > | and 2 — > zy. When fci, fc 2 and 
A; 3 are large, the saddle point for this integrals are x c = y c = z c = 1 ; 
however the denominator vanishes at this point. Therefore, we first deform 
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the contours. The integral of (76) is invariant, except for a sign, under the 
change of variables, a; — > -,u — > - and z — > -, which transforms the contour 

o ? x ' ,y y z' 

around x = y = z = into a contour atx = y = z = oo. By Cauchy 
theorem, the sum of all residues has to vanish if we include the residues 
at infinity, Thereby we obtain the following identity between the different 
contour integrals for (76). 

F = -\[i + <f +1 ]F (77) 

x=y=z=0 A Jx=l,y=0,z=Q Jz=l,x=0,y=0 Jx=±,y=0,z=0 

where F is the integrand of (76). This identity is derived from the invariance, 
except for the overall sign, under the change x — > \,y — > \-> z ~ > \ f° r the 
expression (76). We have 

dxdydz (x + y ) 2kl ( 1 + y ) 2fc2 ( 1 + zy) 2ks 



x= y =z=o (2-Ki) 3 x k iy k i+ k i +k 3 +1 z k 3(l -x)(l- z)(l - xz) 
dxdydz (x + y ) 2kl ( 1 + y ) 2k2 ( 1 + zy) 2ka 

x=y=z=oo (27T«) 3 ^1^^1+^+^3+1^3(1 - x)(l - Z)(l - XZ) 

(78) 

Therefore, we obtain the identity (77) with the factor 1/2. 

The double pole , which appears for the contour integral around 

x — 1, is transformed into a single pole by an integration by parts over x, 
and the sigularity at z — 1 is cancelled by the numerator. Therefore we may 
now use the saddle point at z c = 1, and obtain the large ki behavior, 



+ yjhifa + h) - yjhih + fa)] (79) 

In order to appreciate the asymptotic behavior of these integrals, we have 
compared the exact value for the contour integral (76), multiplied by a factor 
(2k 1 )\(2k 2 y.(2ks)\ : to the large k estimates. For instance for k\ = 30, k 2 = 
15, k 3 = 50 the exact integral is 7.4632 x 10 55 , and the large ki asymptotic 
formula (79) gives instead 7.4864 x 10 55 . 

We have considered the region u\ > u 2 > u^. There are other regions Ui > 
Uj > ui. Adding their contributions amounts to summing over permutations 
of the ki, k 2 and k 3 ; taking also into account the I kl > k2 > k3 for h odd, we obtain 

Aki+k2+k 3 

- 2^!(2^)!(2 fc 3)!(^^ 2 ^3) [ ^^ + 
+ ^k 2 (h + k 3 )} (80) 
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This leads to, 



X; < tiM 2kl trM 2k2 trM 2ki > c 
N 3 



-[\jh{h + k 2 ) + Jhfo + fa) + yjk 2 (k 1 + k 3 )] 



27r(k l + k 2 + k 3 ) 

(81) 

We have evaluated the leading term of order one in the large N limit of 
the three-point correlation function. However, this leading term is cancelled 
when we consider connected correlation functions. Indeed let us consider the 
expansion of U(t) and U(t\,t 2 ) : 

rr/,N I f du N . it . N _£_ 

ult) = Nt^iu e ^ + m> e " 



it2\N 



^ 

Then 



U t 2 ) = J-e-*(M) f ^e^^{l + ^ni + ^. 
C[ U 2) N* 6 f (2niY (u 1 -u 2 + ^)(u 2 -u 1 + f: 

= ^t 1 - + + ^1*2 + + t 2 f - ^{t\t 2 + #| + tl t\) 

^(^ 2 + ^) + ^2 + 0(t 6 )] (83) 



x 



J 11 (2m) ( Ul - M 2 + f )(« 2 - M 3 + f )(«3 - «1 + f ) 



^2 [1 - +t 2 + hf + ^(tx + t 2 + t 3 ) 4 + -^(tl*2 + tlt 3 + * 2 * 3 ) 

~2^(*?*i + + *!*!) - ^(^ 2 + ^ + + ^ 3 + tit| + 1 2 *|) 

-^t^tsih +t 2 + t 3 ) + ^(t! + t 2 + hf 

+^t 1 t 2 t 3 (t 1 + t 2 + t 3 ) + 0(t 6 )] (84) 

Combining these expansions one obtains 

U c (h +t 2 + t 3 ) - NU (*! + t 2 , t 3 ) - iVf/ (h + t 3 , i 2 ) - NU (t 2 +t 3 ,t ± ) 
+N 2 U c (t u t 2 ,t 3 ) 
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= ^ [^(*? + 4 + 4) + \(44 + + 44) + *i* 2 t 3 (ti + * 2 + * 3 ) 

+ g(*?*2 + t?<3 + ^1 + ^3 + <3*1 + *3*2) + 0(t e )] (85) 

It is order of Thus we see that the term, which we had considered in the 
large N limit, is cancelled by the additional terms terms in (85). In other 
words, the odd-power of are cancelled in the combination of (85). We thus 
have to expand also the denominator in (84), and compute the order , 
instead of the order in U c (t 1 , t 2 , t 3 ). 

Noting that expressions for the large N limit of < trM 2kl tiM 2k2 > is 
given by 

< trM 2fel trM 2fe2 >~ A M \ (86) 

N 2 7rN 2 (k 1 + k 2 ) v 1 



we have 



< trM 2kl+2k2 trM 2k3 > +— < trM 2kl+2ka trM 2k2 > 
N 2 N 2 

+ < trM 2k2+2ks trM 2kl > 
N 2 



h(h + k 2 ) + 



irN 2 (k 1 + k 2 + k 3 )> 

^Jk 2 (h + h) + \/h(k 2 + h)} (87) 

This sum is identical to what we found in (81) for < trM 2kl trM 2k2 trM 2k3 >. 

To discuss the next order terms, we first derive a formula for the cor- 
relation functions of the vertices, which is applicable to the general case 
of n-point vertex correlations for arbitrary genus. Let us use the following 
notation, 

[hk 3 ] = ^ (88) 

Returning to the expression for the correlation function of n-point vertices 
in the large ki limit, one shifts Ui — » u\ — 

V it-k 

tt(+ + \ -^Vfc 3 I Uduj " "i 

= G T^ E k 't I J^ e J2 ( V _thM}^L) . . . ( V HM*]) 1 * ) 

j (2tt<)» («i - ^) ui+1 tr k - «i)^ +i ) 

(89) 

From this representation one extracts < Il"=i trM 2 ^ > c as coefficient of 
the relevant power of tj. 
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The above formula is applicable to the n-point case for arbitrary genus. 
For example, we consider n=3, and evaluate < llf=i tiM 2ki >. 
We use the notation, 



with C/(0) = 1. Expanding 



1 ^ifc^y (9i) 



the integral becomes 



/ M l u 2 u l 

J ai (2t 1 )J n (2t 2 )J 83 (2t 3 ) (92) 

The u-integrals are now easy, and the coefficient of (iti) kl {it 2 ) k2 {it 3 ) 2k3 , de- 
noted I^^ k J 3 ' 2k3 , is expressed as a sum. This sum is given by the contour 
integral 

iPyf/ ! ^E(*»)*(»)»(j)*C ll (^)Q,Wa.W (») 

We consider separately (i) /i,/ 2 even, and Z 3 odd, and (ii)/ 1? / 3 even, and l 2 
odd. When we consider the universal scaling limit for large hi, this difference 
can be neglected. We replace U — > 2/j or U — > 2^ + 1. The sum over /j 
becomes 

E^MW ~EM")W - (1 _ 2 J y) , +1 (94) 

For instance when z/ = 3, we have 



Y,C 2l (3)(xy) 1 = i^(2/ + 3)(2/ + 2)(2/ + l)(^) i 
1 + 6xy + (a;?/) 2 



[1-xyY 



(95) 



and, since the saddle point is x c = y c = 1, the numerator at x = y = 1 is 
indeed equal to 2 3 . 
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Then one has 

C5f ^ = -(~l) V1+Va [klk2] V1 [k2k 3 \ V *[hk 1 ] V3 

• dxdydz (1 + x) 2kl (l + y) 2k2 {l + z) 2fc3 2 Vl+v * +V3 



x 



11(2^)! / (27ri) 3 x ki-^-^-yk 2 +^ +1^3(1 _ x^^+^l - § ) iy2+1 ( 1 - a*) 1 ** 1 

1 r dxdydz (-l) U3 (x + yf kl {l + y) 2k2 (l + yz ) 2k3 2 Ul+U2+u * 
' Il(2fci)! * (27ri) 3 a;fciyfci+fe+fe3+i^3( x _ 1)^+1(2 - - 1)^+1 

xtfcifel^tfcafcsl^tfcafci]^ (96) 

In the last line, we have changed x — » | and z — > zy. We have also dropped 
the subleading powers of re and y for large k{. 
We have from (??), 

-L < ntrM 2 ^ > c = n(2%)!(E«f fc3 )^ E ^^ (97) 

As discussed before, we need the term of order of -^4, and the order A; 3 / 2 
in the large k limit. For this reason, we take v 1 + u 2 + u 3 = 1 in (96). 
We have 



1 1,1 1 



(z — l)(xz — 1) z(x — 1) 2 — 1 xz — V 
Using this identity, we have for v\ + v 2 + vj, — 1, 



(98) 



AT (x-l) 2 (2-l)(x^-l) (x- l)(z- l) 2 (xz- V 

ihh] 

(x- \){z-l){xz- 1) 



2 



A; 3 k 2 + ks k-i k 3 + k\ , . 



iVs?s 3 2iVs 2 s 2 JVsf(si + s 3 ) 2iVs 2 ( Sl + S3) 1 

where we have expanded x, y and z near the saddle points as x = 1 + isi, y = 
1 + 1S2, z = 1 + is 3 . From the saddle point analysis, we obtain 

Iff] ^ e -^4- tl± ^4- ! f4+^+ 1 f^sj 

NJf = \(2m) e 

= h + h (100) 

where l\ corresponds to the first two terms of J, and I2 corresponds to the 
third and fourth. For the s-integral, we integrate by parts, which reduces the 
integrand of J to the sum of constant terms and single pole terms in si and 
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S3- We take only the constant terms after diagonalization of the quadratic 
form in the exponent. We obtain 



y/kik 2 k 3 (k 2 + h) 

1 ~ 8A^7T 3 /2(A; 1 + k 2 + fc 3 ) 

t = Vhhhih + fc 3 ) , . 

2 SAvr 3 / 2 ^ + A; 2 + fc 3 ) 1 J 

Returning to the expansion of (92), we add the permutations over the k iy 
and obtain 



ij < tr*»trJK«trM* > c = >™ 
Dividing by by 4kl+k2+ ^Vkik 2 k 3 ^ we 0D ^ a j n ^ e intersection number 



< > 9 =o= 1 (103) 

The intersection numbers, for the three point n — 3, for higher genuses 
are evaluated by considering higher values for z/j in (97). After factoring 

out 4 fc i+*2+*3Vfcifc 2 fc a> the vertex correlations ^3 < trM 2fel trM 2fe2 trM 2fca > c 

become polynomials in and scale as j^. The intersection numbers are 
obtained from the coefficients of these polynomials. 



5 Characteristic polynomials for Airy matrix 
functions and the replica method 

We consider in this section the correlation functions for the characteristic 
polynomials of random matrices. The Airy matrix model of Kontsevich 
type, will then be derived from these correlation functions at the edge of 
the spectrum. 

In a recent work we have studied the average of products of characteristic 
polynomials [7], defined as 

k 

F fc (Ai,...,A fc ) = < Y\_ det(A a — M) >a,m 

a=l 




dM J] det(Ai • I - M)e"f trM2 + NtrMA (104) 



where M is an x N Hermitian random matrix. It was shown that this 
correlation function has also a dual expression. This duality interchanges N, 
the size of the random matrix, with k, the number of points in F k , as well as 
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the matrix source A with the diagonal matrix A = diag(Ai, A&). Indeed 
we had derived that this same correlation function is given by [7] 

F k (X 1 , X k ) = J dSn[det(a j -i J B)]e-f trB2+iNtrBA (105) 

where A = diag(Ai, X k ) and B is a k x k Hermitian matrix. 

If we specialize this formula to a source A equal to the unit matrix, 
providing thus a trivial constant shift for M, the formula (105) involves 

det(l -iB) N = exp[Mrln(l - iB)] 

= exp[-iNtrB + ^-trB 2 + i^-trB 3 + ■ ■ ■] (106) 

The linear term in B in (106), combined with the linear term of the 
exponent of (105), shifts A by one. The B 2 terms in (105) cancel. In a 
scale in which the initials A^ are close to one, or more precisely N 2 ^(X k — 1) 
is finite, the large N asymptotics of (105) is given by matrices B of order 
jy-i/3_ Then the higher terms in (106) are negligible and we are left with 
terms linear and cubic in the exponent, namely 

r N 

F k (\ u A fc ) = / dB J] e *>* 3 ™^-D. (107) 

This is nearly identical to the matrix Airy integral, namely Kontsevich's 
model [15], which gives the intersection numbers of moduli of curves. 
The original Kontsevich partition function was defined as 



Z 



= ±_J rfMe -itrAM2 + i tr M3 (10g) 



where Z' = J dMe~2 trM . The shift M ■ M - / A. eliminates the M 2 term 
and one recovers (107). 

Let us examine the simple case of a 2 x 2 matrix M, with two points Ai 
and A 2 . (Indeed this simple N = 2 case is useful as a check for the intersection 
numbers. For higher N one could perform a similar analysis.) Then we have 

Z = v /a^(A 1 + X 2 )e^ +X ^Y (109) 

where 

Y = JdMe^ A2M+i ^ M3 

= f dxidx2 ^lf^l )e i (A?x 1+ A|, 2)+ i(,3 +a; 3) 

J A 1 — A 2 

)[^=e-^ A ? +A ^^(A 1 )^(A 2 )] (110) 



A 2 

d d 



Xl-X^dXl dX 2 2 n ^X 
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where 



/ dxe x ° x+ f- xi 



Then we get 



*(A) = - g- (HI) 
J e 2 A dx 

Z = 1 + l i3 o + k ii+0{ h ) > (112) 

log Z = ^ + ^ + ^ + 1*? + lf f 2 

1 2 ^ 2 1 ~3 1 ~ ^ 1 ~ ~ ~ 

+ gMo + y2*i + 43^^ + j2 totlt2 

+ ik h+0< ^ (113) 

where we have defined the moduli parameter U as 

11-11- 1 1 

*0 - Y" + T' 1 ~~ \3 + T3' tn ~ T^+I + T2^+T- ( 114 ) 
Ai A 2 A x A 2 A 2 

(We need the notation t n , instead of the usual notation t n , for distinguishing 
those parameters from the Fourier transform parameters U). 
From these coefficients, we recover the known results 

< Tq > ff =o= 1, <r 1 > g=1 = — , <r V 1 > 9=0 =l, 

2 1 1 

<T l > g=1 =— , <r r 2 > 9=1 = — 

1 1 

< r o T i ><?=o= 2, < rj 3 > s= i= — , < r 2 r 3 > 9= i= — 

< T TiT 2 > # = 1 = < r 4 > s=2 = — L (115) 
The generating function for the intersection numbers < II > 3 is 

oo Tm/ 

bgz = E<Cn m2 >n— r ( 116 ) 

mi 1=0 mi - 

where is related to the eigenvalues of the matrix A as 

^E^# (117) 

Note that the generating function of the intersection numbers (116) is 
quite different from the generating function (67), which was obtained from 
the correlation functions of the vertices. The relation between ti and Xj are 
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where Xj = 7^73, in which kj is the power in M 2fcj . However the degenerate 
case needs some precaution, when we have tf ~ x l jX l m for j 7^ m. 

The replica method for the correlation functions has been used earlier in 
random matrix theory for the GUE [8, 9]. Following this replica analysis we 
study now two types of correlation functions. The first one is the correlation 
function for the eigenvalues of the matrix M : 

k -y 

p(Ai, A fc ) =< n -jrM{K - M) > a,m (119) 

The second type of correlation functions is the average of products of char- 
acteristic polynomials [7], 

k 

F fc (Ai, A fc ) =< 1] det ( A « " M ) >a,m ■ (120) 

a=l 

where the average < ... >am is with respect to the probability distribution 
Pa{M), 

p A ( M ) = L e -^rM* + NtrMA (m) 

z 

The random matrix M is a complex Hermitian N by N matrix and A is an 
external source, which we can take as a diagonal matrix, A = diag(a 1; ajv) 
since the integration measure is unitary invariant. If A is zero, it reduces 
to the Gaussian unitary ensemble (GUE), but it is convenient to use the 
probability distribution Pa{M) for setting up the replica method, even if we 
let A = at the end. The correlation functions with the distribution Pa(M) 
require the HarishChandra-Itzykson-Zuber formula [11, 12] for the unitary 
group integral. 

We have evaluated the F.T. of the correlation functions near the edge in 
the previous sections. To prove that F.T. of the correlation function is the 
generating function of the intersection numbers, we express these functions 
as the zero replica limit of characteristic polynomials average. Let us begin 
with the one point function, namely the density of states ; we first use the 
identity 

X X c) 

< tr5(\ - M) > a ,m= 1™ < t det (^ - it ~ M )T >a,m (122) 

and use the duality derived in [7] to write 

N 

< [det(A - M)] n > a ,m=< II t det K - iB)] >a,b (123) 

7=1 

where B is an n x n random Hermitian matrix, and A, in this case, is a 
multiple of the n x n identity matrix : A = diag(A, A). Note that we have 
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traded an n-point function of N x N matrices M for an N point function of 
n x n matrices B. 

If we explore the edge of the distribution by taking again a 7 = 1, and 
expand in powers of B, we find that the B 2 term cancels, and obtain 

U(t) =< —tre ltM >= 3m lim / d\e itx --^- [ ^g^f b 3 +*b(a-i)] (124) 

After integration by parts over A we obtain 

U(t) = t lim| dAe itA i J dBe itr lTB 3 +iB(A-i)] (125) 

Since the replica parameter n means the repetition n times of the same A , 
we simply replace the t\ of the previous un-replicated case by 

k 1 k n 

j=l A i j=l A j 

For the one point function k=l, the expansion of the Airy matrix model 
in terms of t\ is thus a power series in n. In the zero-replica limit one can 
neglect all terms beyond the linear one in n. Thus, U(t) in (125) is expressed 
as linear combinations of the ti. The F.T. of ti — tt yields dimensionally a 
factor t l . Therefore, U(t) appears as a power series in t, whose coefficients 
are the intersection numbers < 77 >, the same coefficients obtained in the 
Kontsevich model. 

This argument holds also for the two point correlation function, since 

lim 7— -—7— — -jr < [det(Ai - M)det(A 2 - M)] n >= tr - 1 - tr - 1 - (127) 
n^o 9Ai9A 2 n 2 1 K ' K n Ai - M A 2 - M 

In this two-point case (k=2), one deals with Ai and A 2 and i\ becomes 

tl = A[ + A| (128) 

The replica limit (n — > 0) requires to retain the terms of order n 2 , i.e. 
products of the form t\t m . Then ?7(ti,t 2 ) is a power series in ti,t 2 , whose 
coefficients are the intersection numbers < rir m >. 

The F.T. of the k-point correlation function U(t\, ...,tfc) produces, by the 
same replica method, the intersection numbers < r mi r m2 ... > as coefficients 

ofc^r-'-c- 

Thus we have shown that the F.T. of the correlation functions in the 
edge region, is the generating function of the intersection numbers , our 
main conclusion. 
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6 Summary 



We have derived the expressions for the correlation functions of the vertices 
< IliLi trM 2fei >, when N and fcj are simultaneously large, in a scaling 
region in which k 3 scales like N 2 . The expressions are obtained from an 
exact contour integral representation valid for finite N and for an arbitrary 
external matrix source. The coefficients of the expansion of the correlation 
functions provide the intersection numbers of the moduli space of curves. 

The correlation functions of the eigenvalues may also be obtained from 
the average of characteristic polynomials in a zero- replica limit. Using a 
previously derived duality in which the size N of the matrix is interchanged 
with the number of points of the correlation functions, we have recovered the 
Airy matrix model of Kontsevich and re-derived the intersection numbers by 
a simple saddle-point analysis. 
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